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Abstract: We have previously described a mathematical formulation for a parton shower 
based on the approximation of strongly ordered virtualities of successive parton splittings. 
Quantum interference, including interference among different color and spin states, was 
included. A practical numerical implementation strategy was left unspecified. In a sub- 
sequent paper, we showed that if we add the further approximations of taking only the 
leading color limit and averaging over spins, we obtain a shower evolution that can be 
implemented as a Markov process. In this paper, we outline a strategy for including the 
correlations induced by parton spins. 
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1. Introduction 

In Ref. [1], we described how the evolution of a leading order parton shower could be 
formulated so that the shower is built on the approximation that the virtuality in each 
successive parton splitting is much smaller than the virtuality in the preceding splitting. 
No further approximation is made. The evolution equations thus necessarily account for 
interference graphs, for the color state of the partons, and for the spin carried by each 
parton. The evolution equations generate a desired cross section as a nested set of integrals 
that could, in principle, be performed by numerical integration. However, we did not 
present a method for implementing the integrations in a practical manner. 

Within the general framework of Ref. [1], it is possible to make further approximations. 
In particular, one can take the leading color limit, l/A'^ — > 0, where A'^c = 3 is the number 
of colors. This greatly simplifies the structure of the evolution. Furthermore, at each 
splitting one can average over the spin of the mother parton and sum over the spins of 
the daughter partons. Both of these approximations are commonly used in parton shower 
Monte Carlo event generators. Here, we still keep interference graphs, as in those event 
generators that are based on color dipoles. In Ref. [2], we saw that the general formalism 
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with these further approximations can be formulated in a standard style of calculation as 
a Markov process. 

What might be the effect of averaging over spins? Partons carry quantum spins and 
evolution creates multi-parton states in which the spins of the partons are correlated with 
one another. At each successive splitting, the spin of the mother parton can affect the 
distribution of the momenta of the daughter partons. In this paper, we assume that the 
spins of the final partons that appear at the end of the shower are not measured. That 
is, we take the square of the amplitude to produce each spin state and we sum this over 
the spins. Nevertheless, the spin distribution of the intermediate partons is imprinted on 
the momentum distribution of the partons. Thus, to get the momentum distribution right, 
one needs to follow the spins. 

Additionally, parton spin correlations are important for particles that decay via their 
coupling to VF-bosons and, in general, for beyond-the-standard-model particles, as de- 
scribed in Ref. [3]. For this reason also, one wants to include spin effects in a QCD parton 
shower in a simple way that can be extended to other interactions. Finally, perturbative 
QCD matrix elements contain the full spin information for the intermediate virtual parti- 
cles. If one wants a shower to match these matrix elements in the soft and collinear limits, 
then the shower should also contain the full spin information for the intermediate partons. 

The purpose of this paper is to study how, still working in the leading color approx- 
imation, one can put back the effect of the parton spin states on the final momentum 
distribution. Suppose that one starts with a spin averaged shower. Then there is a weight 
factor that should be associated with each event that gives the probability to get this event 
with spin included divided by the probability to get this event (with the same splitting 
history) in a spin averaged shower. Our approach will be to calculate the weight factor 
and associate it with the event, so that the cross section to produce a given configuration 
of partons is proportional to the sum of the weights of the events for which the partons are 
in that configuration. 

It might seem that a calculation of the weight factor is very complicated, involving, 
as it does, the entangled spins of states with many partons. However, following an insight 
in a work by J. Collins [4], we find that the calculation of the spin weight factor is quite 
straightforward and uses an amount of computational resources, both memory and time, 
that is proportional to the number of partons in the event. 

The approach of this paper can be compared to that of the two most commonly used 
parton shower event generators. In Pythia [5], one simply averages over spins. However, 
some spin induced correlations are incorporated by letting the azimuthal angle of each 
splitting be correlated with azimuthal angle of the splitting that produced the mother par- 
ton, following the prescription of Ref. [6] . In Herwig [7] , spin is included in a fashion that 
produces the spin correlations in the limit that the splittings are collinear or soft x collinear. 
This leaves out purely soft splittings, for which the angular ordering prescription of Her- 
wig is approximate with respect to azimuthal angle correlations even after averaging over 
spin. The method [4, 8, 3] has the advantage of not requiring weight functions. However, 
the method is not consistent with the use of exact momentum conservation at each splitting 
nor is it compatible with the inclusion of the exact angular correlations arising from soft 
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gluon interference diagrams. Each of these is important in the formulation of Refs. [1, 2]. 

In the fohowing section, we describe how spin appears if we view it as adjoined to a 
spin averaged shower. This leads to the definition of the spin weight factor. In Sec. 3, which 
forms the heart of this paper, we describe how to evaluate the spin weight factor. Then in 
Sec. 4 we investigate some properties of this factor. In Sec. 5 we give a numerical example 
of how one step in the evolution of the spin weight factor works and in Sec. 6 we give a 
numerical example of how three steps in this evolution work together to correlate azimuthal 
angles of parton splittings. In Sec. 7, we turn to a technical topic, how we partition the 
coherent radiation from two partons, call them I and k, into two terms, one treated as 
radiation from parton / and one treated from parton k. Our treatment of the partitioning 
function Ai^ generalizes our earlier work on this function, allowing, in particular, an An. 
that is spin dependent. We present several possible choices for A/^, any of which are 
compatible with the methods of sections 2 through 6. The numerical examples of sections 
5 and 6 are based on one of these choices. Finally, Sec. 8 contains some concluding remarks. 

2. Incorporating spin in the parton shower 

In this section we take the general parton shower formalism of Ref. [1] and make the 
leading color approximation at each step of the shower but retain the spin information. In 
Ref. [2], we have seen how to formulate a leading color, spin averaged shower. We write 
the operators that occur in the leading color shower with spin as products of the functions 
that apply for the spin-averaged shower with specific operators that act on the partonic 
spin space. This will leave us with a matrix element in the partonic spin space. This allows 
us to express the parton shower with spin using the parton shower without spin and then 
incorporating spin as a weight factor. In the following section, we will see how this spin 
matrix element can be evaluated. 

2.1 Parton splitting with spin 

We consider a shower in the leading color approximation, but keeping the quantum spins. 
We start with the complete shower evolution as described in Ref. [1], but we make the 
leading color approximation as described in Ref. [2] . We need states representing a set of 
m final state partons and two initial state partons with momenta p and flavors / in a color 
state c specified by a color string configuration. Our states are further specified by two 
spin indices, s and s', for each parton. The spin index s describes the spin of that parton 
in the quantum amplitude |M) while the spin index s' describes the spin of the parton in 
the conjugate amplitude {M\. We call this state \{p, f,c, s' , s}m) ■ Our notation is that 
the partons carry labels a, b, 1, . . . , m and that, for instance, {p}m denotes the ordered set 
of momenta {pa,Ph,Pi, ■ ■ ■ ,Pm}- 

The states evolve with a linear operator U^'^{t{,t') that acts on the space for which 
the \{p, f, c, s', s}m) are basis states. Here the superscript "Ic" designates quantities in the 
leading color approximation. The evolution parameter is a shower time t, the logarithm 
of the virtuality in parton splittings. The operator L{^^{ti,t') effects the evolution from 
a time t' to a final time t{ at which perturbative shower development is halted. The 
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evolution equation for lA^'^ can be specified by stating the action of lA^'^ on a general state 
/, c, s', with m — 1 final state partons/ 

U''{ti,t')\{p,f,c,s',sU-i) = 

A(°) (tf , t'; {jl, /, c}m-i) I {p, f, c, s', s}m_i) 

+ t'dr^^ I [d{pJrcU] U''{k,r)\{p,f,c,s',s}m) ^^'^^ 

X {{p, f, c, s', sUlnHr) I {p, /, c, s', sU^i) A(0)(t, t'; {p, /, c}^„i) . 

The factor A^^\T,t';{p,f,c}m) is the Sudakov factor that gives the probability for the 
state not to undergo a parton splitting between shower times t' and r. The Sudakov factor 
is the same as in the leading color, spin averaged approximation, which we designate with 
a superscript (0). Thus the first term in Eq. (2.1) represents the possibility that the system 
evolves from t to tf with no splitting. In the next term, the state evolves from t to r with 
no splitting, then undergoes a splitting according to the splitting operator H^f. After the 
splitting there are m final state partons, with the new parton carrying the label m. We 
integrate over the time r. We also integrate and sum over the momenta, flavors, and colors 
of the new partons after the splitting, as in the shower averaged over spin. Since we now 
include spin, we also need to sum over the spin indices {s',s}m of the partons after the 
splitting. 

Using Refs. [1] and [2], we can write the matrix elements of the splitting operator in 
the leading color approximation as 

({P, /, C, S, s}rn\'H^l{t)\{p, f, C, s' , s}rn-l) 

nc{a)nc{b) r]a,ri\, /a/yl(^a, lJ^F)fb/B^'^^-> /^f) 



Tfl 



E 

X {{P, fU\ri\{p, f}m-l) S{t - Ti{{p, /}, 

x{0{fm = g) Y {{^'^^}rn\y{l,k-,{f,pU)\{s',sU_i)'^ik{{p,f}r 
V I, 



k 



X {{c}m.\ali^\{c}m~l) 
+ 0{L^g) {{s,sU\y{l,l;{f,pU)\{s',sU) <^ll{{pj}m) 

X e(/zG{g,g})({c}™|at(/)|{c}„_i)+e(/i = g)({c}^|at(/)|{c} 



m— 1 / 

■(2.2) 

The first line on the right hand side of this formula contains factors copied directly from 
Ref. [1]. There is a sum over the index / of the parton that splits. Then there is a ratio 
of parton distribution functions, momentum fractions 77, and the numbers of colors ric 



^In Refs. [1] and [2], we started with a state of m partons, so that the spUtting resulted in a state with 
m + 1 final state partons. In this paper it is most convenient to decrease m by one. 
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carried by the partons. This ratio is 1 for a final state splitting but different from 1 for 
an initial state splitting.^ The next line concerns the relation of the variables {p, f}m and 
t to the variables {p,f}m-i- For the flavors, this factor vanishes unless fm + fi = fu 
with the evident definition of adding flavors, and it vanishes unless fj = fj for the other 
partons. For an allowed relationship between {/}m and {f}m-i, the flavor factor is 1. 
There is a similar factor for the momenta. Given the momenta {p}m-i, the momenta {p}m 
must lie on a certain three dimensional surface specified by the momentum mapping TZi 
defined in Ref. [1]. The function {{p, f}m\Vi\{p, f}m-i) contains a delta function on this 
surface. There is also a delta function that defines the shower time t as the logarithm of 
the virtuality in the splitting, 

Ti{{pJU)=log(- ^ ,^ . ] . (2.3) 

\\{pi + i-iy^'^+^^'^Pm)^ - m?{fi + fm)\J 

Thus if we integrate [{p, f,c, s' , s}m\'H^{t)\{p, f,c, s' , s}m-i) over t and the momenta 
{p}m, we are really integrating over three variables that describe the splitting of parton /. 

The most important functions in Eq. (2.2) are the functions y{l,k;{f,p}m), which 
describe parton splitting with gluon emission. These have the structure 

y{l,k;{f,pU) = {w(/,/;{/,pW 

^lk[{pj]m) I (^2.4) 
-Alkmm) \w{l,k-{f ,p]r,^) + W{k,l-{f ,p]rr. ' 



The functions W(/, Z; {/,p}m), yV(/, A;; {/,p}m)5 and Vy(A;, /; {/,p}m) are operators on the 
combined parton spin space. Thus we can specify them by giving their matrix elements 

({s', s}„^|>V(/, k; {f,p}m)\{s\ S}m) ■ 

These operators are fully defined in Ref. [1] and there is little point in repeating the 
definitions here. However, the physical meaning can be appreciated by using some simple 
pictures. 

The spin matrix element of W(/, /; {/,^j}m) is illustrated in Fig. 1. In the process 
described by this function, one of the partons, with label splits in the quantum amplitude 
to form partons with labels I and m. In the complex conjugate amplitude, parton I also 
splits to form partons / and m. For each parton with label i in the amplitude, the spin 
before splitting is Si and the spin after splitting is Sj. For each parton with label i in the 
conjugate amplitude, the spin before splitting is s'^ and the spin after splitting is s^. The 
function W(/, {/,p}m) has a nontrivial dependence on the spins of the active partons I 
and m, a dependence that is taken from the Feynman rules for the splitting amplitudes. 
There are several other partons, "1," "2," etc, that are simply spectators. For a spectator 
parton, W(Z, I] {/,p}m) contains factors ^s^^sj and 5^' that keep the before and after spins 
the same. 



^Initial state splittings are done with backwards evolution. Our notation is that the momenta pa and pb 
of the initial state partons denote their physical momenta, whereas the flavors /a and /b denote the flavors 
leaving the hard interaction, which are the opposite of the physical flavors entering the hard interaction. 
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Figure 1: Illustration of the function Wu in Eq. (2.4). Parton m is emitted from parton I in the 
quantum amplitude and from parton / in the conjugate amplitude. 

The spin matrix elements of W{1, k; {f ,p}m) and W{k,l;{f,p}m) are illustrated in 
Fig. 2. These functions describe quantum interference between two graphs. Consider 
the function k; {f ,p}m)- In the quantum amplitude, parton / splits to form partons 
with labels / and m. In the complex conjugate amplitude, parton k splits to form par- 
tons k and m. There are several other partons that are simply spectators. The function 
W{l,k;{f,p}m) is proportional to the unit operator on the spins of partons / and k but 
has a nontrivial dependence on the spins of parton m. This spin dependence is taken from 
the Feynman rules for the splitting amplitudes making use of the eikonal approximation. 
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Figure 2: Illustration of the functions Wik and Wki in Eq. (2.4). The left hand diagram illustrates 
Wik, in which parton m is emitted from parton I in the quantum amplitude and from parton k in 
the conjugate amplitude. The right hand diagram illustrates Wki, in which parton m is emitted 
from parton k in the quantum amplitude and from parton I in the conjugate amplitude. 

In Eq. (2.4), the operator Aik tells how the l-k interference graphs are partitioned 
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between a term treated as a splitting of parton /, with parton k playing an auxiliary role, 
and a term treated as a splitting of parton k, with parton / playing an auxiliary role. In this 
paper, we allow Aik to depend on the spin indices of parton m, so that it is an operator on 
the parton spin space. In Ref. [2], we specified Aik{{p}m) as a spin independent function 
of the momenta. We give possibilities for Aik in Sec. 7. The structure of the calculation 
outlined in this and the following section does not depend on which of these choices of Aik 
one takes, although certainly the numerical performance of the algorithm can depend on 
the choice. We when we provide numerical examples in Sees. 5 and 6, we make a definite 
choice that will be explained in Sec. 7. 

In Eq. (2.2), the spin matrix element of y{l, k; {f,p}m.) multiplies a spin independent 
function ^ik- This function appears also in the denominator of y{l,k; {f,p}m), so that it 
cancels. The function is the spin averaged splitting function, the spin average of the 
spin dependent splitting functions in y{l, k; {f,p}m), defined by setting {s'}m = {s}m and 
summing over these after-splitting spins and setting {s'}m-i = {s}m-i and averaging over 
these before-splitting spins. The spin matrix elements of 3^ are the ratios of the splitting 
functions with spin to the spin- aver aged splitting functions. Thus they are normalized to 

E {{^^^}m\y{l,k;{f,pU)\{s,s}m)=l . (2.5) 

{s}m {s}m-l 

The functions are given in Ref. [2]. They can be recovered from the matrix elements 
of y by using the normalization condition (2.5). 

The spin matrix element ofy{l, k; {f,p}m.) multiplies a color factor ({c}m|a|fc|{c}m-i)- 
This factor equals 1 provided two conditions hold. First, partons / and k must be color 
connected in the initial color state {c}m-i- Second, the new color state {c}m must be the 
same as {c}m-i with the gluon with label m inserted between partons / and k. If either of 
these conditions fails, this factor vanishes. 

This completes the brief description of the factors in the term containing y{l, k; {f,p}m) 
in Eq. (2.2). There is a second term, which covers the possibility that the newly created 
final state parton, labeled m, is not a gluon. Then there is no interference diagram to con- 
sider so there is no second parton with label k. Our notation for this case uses functions 
with the same names as for when there is a helper parton k, but sets k ^ I. One can have 
/m 7^ g in two cases. In the first case, one has an initial state splitting in which {/;, fm} 
is either {g, g,g} or {^, g, g}. The color factor for this splitting, ({c}m|og(0 |{c}m-i), is 
1 if the new color state {c\m is the same as {c}m-i with the end of the string that was 
at quark or antiquark I now terminated at quark or antiquark m and the new gluon with 
label / inserted just next to the end of the string. Otherwise, this factor vanishes. In the 
second case, one has an initial or final state splitting in which {/;, /„} is either {g, g, q\ 
or {g,q,q}. The color factor ({c}m,|ag(/) |{c}m_i) is 1 if the color state {c}m is related to 
{c}m-i by cutting the color string on which parton / lies into two strings, terminating at 
the new quark and antiquark. Otherwise, this factor vanishes. 

For k = I, the definition of y{l, I; {f,p}m) is very simple. 



The operators W are specified in Ref. [1]. The color factor C{fi,fm) corresponds to the 
color operators in the two cases in Eq. (2.2) with fm'/'S 

nif f\ {/« Jm} = {g,g} or {g,g} 

[Tr {jl,fm} = {q,q}or{q,q} 

Again, the functions are the spin averaged splitting functions. They are given in Ref. [2]. 
They can be recovered from the matrix elements of 3^ by using the normalization condition 
(2.5) with k = l. 

We can think of the operators 3^ as generating evolution in spin. They operate on 
the space of spin states with two spin indices for each parton, with basis states |{,s', s}m)- 
These operators map the spin state for m — 1 final state partons into a spin state for m 
final state partons. In general, the spin states are linear combinations, |/9spin)i of the basis 
states, 

\&= E \W^^U{{s',sU\p^-l) . (2.8) 

{s',s}m 

We will use subscripts "spin" for vectors in the spin space other than the basis vectors 
s}m) in order to emphasize that these are vectors in the finite dimensional space with 
2x2 dimensions for each parton. 

2.2 Starting point with spin 

The starting point for evolution is a state for two final state partons and two initial state 
partons, assuming that we start with a 2 — > 2 hard process. Such a state is a mixture 
of the basis states \{p, f, c, s', s}2) that represent a complete description of the momenta, 
fiavors, colors and spins for the partons, including two spin indices but, since we work in 
the leading color approximation, only one color index. The leading color starting state is 
then |p''^(t)) at shower time t = and is represented as a linear combination of basis states 
as 

W"iO)) = ^ [[d{pJ,c}2] \{pJ,c,s',s}2){{p,f,c,s',s}2\p'%0)) . (2.9) 

{S',S}2 

Here ({p, /, c, s', s}2|/j(0)) is obtained from the 2^2 matrix element,^ 

({p, f,c,s, sUp (0)) = (2.10) 

X {M{{p,f}2)\{s',c}2){{s,c}2\M{{pJ}2)) . 

From Ref. [2], we recall that the starting point for evolution averaged over spins (and in 
the leading color approximation) is 

{«}2 



^As explained in Ref. [1], we should most properly project out the component of \M{{p, f}2)) that is 
proportional to a color basis state |{c}2) by using a dual basis state ^({c}2|, but in the leading color limit 
there is no distinction between the dual basis states and the ordinary basis states. 
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It is convenient to define an initial vector in the spin space, /, c}2)spm) that depends 

on the momenta, flavors, and colors of the partons by 



{{s',s}2\p{{p, f, C}2)spin) 



{{P, f, c, s',i 




{{P, f, C}2 


|p(0)(0)) 


{M{{p,f}2) 


\{s',ch){{s,ch\M{{p,fh)) 




\{{s,c}2\M{{p,fh))\' 



(2.12) 



With this notation, we can write the initial state with spin as a product of the initial state 
without spin and a factor that contains an initial vector in the spin space, 

{{p,f,C,s',sh\p"'{0)) = {{PJ,C}2\P^'K0)) X {{s',s}2\p{{p,f,c}2)spin) ■ (2.13) 

2.3 Evolution with spin 

The evolution equation with spin is closely related to the evolution equation averaged over 
spin at each step. To obtain the spin averaged evolution equation from Eq. (2.1), we need to 
eliminate the spin indices and the sum over spins. Then in Eq. (2.2) we need to replace the 
matrix elements of the operators y on the spin space by 1 , leaving only the spin averaged 
splitting functions We can get the full spin dependence back in each evolution step by 
inserting a factor 

Yl \{^''^}m){{S',s}m\y{l,k;{f,p}m) ■ (2.14) 

{s',s}m 

Here the momenta and flavors {f,p}m are the momenta and flavors that appear in the 
integrations for the spinless evolution. In an implementation of the integrations as a Markov 
chain, they are the momenta and flavors chosen at that step. Similarly, I is the index 
designating the parton that split at that step. If the splitting involved the emission of a 
gluon into the final state, then a partner parton with index k is also selected. If not, then 
in that step k = I. 

We have displayed the spin indices in Eq. (2.14), but we can recognize that the com- 
pleteness relation for our basis states allows us to replace 

\{srsU){{s\s}m\ = l . (2.15) 

{s',s}m 

Thus we obtain the full spin dependence by inserting the operator y{l,k]{f^p}rn) corre- 
sponding to the parameters of that splitting. After several steps, the spin state is 

kSi) = 3^(^-> {/,riir^) • • • yilA. h; ) ^^(^3, h; )|p({J>, /, c}2)spin) . 

(2.16) 

Here \p{{p, f,c} 2) spin) is the starting spin state obtained from the hard matrix element. 

(s) 

The first splitting produces 3 final state partons, with momenta and flavors {f,p}\^ , by 
splitting parton ^3 with the participation of partner parton ^3. The number of final state 
partons increases as the shower progresses, so that at the later stage there are m partons. 



- 9 - 



2.4 End of the shower with spin 

At the end of the shower, at shower time t{, we apply an evolution operator that turns the 
partons into hadrons, then measure the hadronic final state with a measurement function 
Fh- Following Ref. [1], we write this as 

^''[^h]= Em / [^{^''^'"}^] {Ft\U''''\oC,t,)\{p,f,c}N) 

^ ' (2.17) 

X {{P,f,C,S,s}N\p^''{ti)) . 

{«}iV 

Here we have made the assumption that neither the hadronization process nor the ulti- 
mate hadronic measurement depends on parton spins. Thus we have computed the total 
probability to get momenta, flavors, and colors {p,f,c}]\f by setting s = s' and summing 
over s. 

At each step of the evolution. 



{{pJ,C,s',sU\p''{tm)) = ({p,/,c}™|pW(t^)) X ({5',s}m|pSn) • (2-18) 

Thus, 



N ^^--^ (2.19) 



X 



^ispinlP'spm; ' 

where 

(lspin|{s',s}m) = n ' (2-20) 

iG{a,b,l,...,r?i} 

SO that 

(lspin|0= (Win|{s',4iv)({.',.}N|p£i)=E({«'^}^lC) • (2.21) 

{s',s}jv {s}n 

According to Eq. (2.19), the leading color cross section including the effects of spin on 
the distribution of partons is the same as the leading color cross section without spin times 
a factor (lspin|/Os^n) • We propose to include this factor as a weight for each Monte Carlo 
event. In the next section, we investigate how to compute it. 

3. Evaluating the spin factor 

We need to evaluate the spin factor in Eq. (2.19), 

(lspin|£i) = (lspin|3^(^JV, kN] {/,P}!?^) • • • 



X y{U, ) 3^(^3, ks; {f,p}f)\p{{p, /, C}2)spin) . 



(3.1) 



This factor contains matrix multiplications, which can be made manifest by inserting spin 
sums using Eq. (2.15). How can we evaluate this spin factor? A little reflection suggests 
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that a Monte Carlo summation, choosing the spin indices at random, is not a promising 
approach. This observation suggests performing ah of the sums exactly. However, j/Og^n) 
is a vector in a space with 2-^(^+^) dimensions. If the number of final state partons is, say, 
100, this vector has more components than any available computer can hold. Thus one 
should be careful about how to arrange the calculation. 

We will follow a method inspired by a work by J. Collins [4] on the subject of including 
spin in parton shower Monte Carlo programs. (See also Refs. [8] and [3].) The problem 
considered in that paper did not include the interference graphs, which we want to include, 
and does not fit well with the way that we have organized a parton shower. However, we 
can use the main idea: that the matrix product in Eq. (3.1) is simplest if we evaluate it 
from left to right. 

To proceed, we denote 

{lspin\yilN,k,,;{f,p}^j^^)...yilm,km;{f,p}t^) = (^ST'^I ' (3-2) 

Then our weight factor is 

(Ispinlpgi) = i^sXlPiiP, /> C}2)spin) • (3.3) 

Following Collins, we can call (^gp™ the decay matrix in spin space. We can generate 
(^plnl recursively, using 

{y}^-'^\ = {yS\yiim, k^-, ifMt^) . (3.4) 

This seems horribly complicated, but it is not. We will see that (^|^^| has the structure 

{yS!\W,^}m)=Uy^riHs'^,S,) , (3.5) 

j 

as illustrated in Fig. 3. That is, the complete decay matrix in spin space is a product of 
decay matrices for the individual partons. Even though (5^^^| has 2^'^"'"'"^) components, 
these components are determined by m + 2 vectors with 4 components each. 

To see why Eq. (3.5) holds, start with the initial condition for the recursion (3.4), 

(nS|{s',4iv) ^ (lspi„|{^',s}^) = n^^;,^. • (3-6) 

j 

This has the promised form^ with 

yg)(4,5,)=5,,,,^. . (3.7) 

Next, we need to see whether this form is maintained under the recursion (3.4). For 
this we need the structure of the spin operators y as given in Eqs. (2.4) and (2.6). This, 

^Eq. (3.7) is the simplest possibility, but one could use a hadronization model that results in a non-trivial 
spin decay matrix for each parton at the end of the shower. For instance, in a string model the polarization 
of a parton at the end of the shower could be correlated with the spatial directions of the string segments 
that couple to the parton. 
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Figure 3: Structure of (i^Jp^^l as given in Eq. (3.5). 

in turn, requires us to examine the structure of the spin operators W as given in Ref. [1] 
and illustrated in Figs. 1 and 2. We first note that the spins of the partons other than I 
and k, that is the partons not involved in the splitting at step m, are left unchanged under 
the splitting: 

yliUi^'p^j) = y^ri\s'^,s,) J i {i^M . (3.8) 

Next, we need to consider the spins of the partons that are involved in the splitting. 

Consider first the case that km = ^m-, which arises when f^a 7^ g- Then we use Eq. (2.6) 
and use the spin operator yV(Z, /; {/,p}m) from Ref. [1]. Writing simply / for 1^.^ the result 
is 



sti: (3.9) 



i' 



X Vi si,si)vl 

{{Pi f}rm Si, Si 

There is a color factor C{fi,fm), Eq. (2.7), and a statistical factor Si{{f}m) equal to 1/2 
for a final state g — > g + g splitting^ and 1 for any other allowed splitting, as given in 
Refs. [1, 2]. In the denominator, there is the spin-averaged splitting function described 
in Ref. [2] . The logic of this is simple. We have the amplitude for a parton of spin si to split 
into partons with spins si and Sm times the complex conjugate amplitude for a parton of 
spin s'l to split into partons with spins s'l and s^. These multiply the spin decay matrices 
for the daughter partons / and m. This equation is illustrated in Fig. 4. 

The case that km / Im, which arises when fm = g, is a little more complicated. We use 
Eq. (2.4) and use the spin operators W{IJ; {f,p}m), W(/,fe; {f,p}m), and W{k,l; {f,p}m) 
from Ref. [1]. To keep the notation simple, we write / for Im and k for km- First, we note 
that the spin of the helper parton is not affected: 



L'li(4,Sfc) = £'^(4,Sfc) • (3.10) 



^The g ^ g + g splitting is part Eq. (3.11) below. 
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Figure 4: Illustration of Eq. (3.9) for spin evolution in a g ^ g + g splitting. 



Second, the new spin matrix for the parton, I, that splits is 



CFSl{{f}m) V- . )\^v(^)(s' sA 

'^ik{{p, /}m) .rt 



Sl,Sl)v1{{p,f}m,s'^,s'i,s[) 

+ 9{le{l,...,m-l},fi = U = g) 

X V21 

- V3,li{p, f}m, Sm, Si, Si) V^i{{p, f}ra, s'm, s'l, S/) 

4:TrasAik{{p}m,s'm^Syn) . r 

^ Osis, Os'S' 

Pm-Pl Pm-Pk ' ' 

;Q)*-Pi e{pm-,s'„^;Q)-Pk 
;Q)*-Pk £{pm,s'^;Q)-pi 



(3.11) 



As in Eq. (3.9), the new spin matrix for parton / is constructed by a matrix multiplication 
from the old spin decay matrix of the daughter partons I and m. There is a color factor Cp 
and a statistical factor Si{{f}m)- In the denominator, there is the spin-averaged splitting 
function described in Ref. [2] . All of this multiplies a factor in braces that is constructed 
from the splitting functions for the quantum amplitudes. The first term is the splitting 
amplitude for a parton of spin si to split into partons with spins si and Sm times the 
complex conjugate amplitude for a parton of spin to split into partons with spins s'l 
and s^. For a final state splitting in which both daughter partons are gluons, there is a 
correction term involving pieces V2 and f 3 of the g ^ g + g amplitude. Following Ref. [1] , 
we write the ggg vertex as the sum of three terms, 

Vi^'^{Pa,Pb,Pc) + V2'^^{pa,Pb,Pc) + Vs^'^ {Pa, Pb, Pc) , (3.12) 



v'^'^'^{Pa,Pb,Pc 



where 



Vi'^'^{Pa,Pb,Pc) = g'^^iPa -Pb)^ 
V2'^'^'iPa,Pb,Pc) = g^^{Pb - Pc)"" 
v'^'^'^{Pa,Pb,Pc) = g^°'{Pc-Pa)'^ 



(3.13) 
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Then vj^i{{p, f}m, Sm, si, si) is the part of the sphtting amphtude built from term J in 
the ggg vertex, vj. This construction from Ref. [1] breaks the symmetry between gluons 
/ and m and ensures that there is a singularity when daughter gluon m is soft, but not 
when daughter gluon / is soft. An alternative is to omit the additional terms and multiply 
both the spin dependent splitting function and the spin averaged splitting function by 
9{z < 1/2). (Here z is the momentum fraction of gluon m, defined, for instance, as in 
Ref. [2].) The remainder of Eq. (3.11) gives the contribution from the l-k interference 
graphs. It is built from the amplitude for soft gluon emission in the eikonal approximation. 
In the eikonal approximation, the spin of parton I remains undisturbed. The function 
^ik{{p}m, s'm, Sm) Specifies how we partition the l-k interference graphs into a fraction Aik 
associated with parton / and a fraction A^i associated with parton k. In this paper, we 
allow Aki to depend on the spin indices Sm and s'^. One, spin independent, choice for this 
function was specified in Ref. [2]. We will give alternative definitions in Sec. 7. The general 
structure of the recursion relation Eq. (3.11) is illustrated in Fig. 5. This structure applies 
independently of the choice of Ai^. 




7 



Figure 5: Illustration of Eq. (3.11) for spin evolution in a. q ^ q + g splitting. The combination 
of the possibly spin dependent An^., the momenta pi and pk, and the gluon polarization vectors for 
gluon m are indicated in the drawing simply as an oval labeled Aik- 

We see that the recursion relation (3.4) preserves the simple structure (3.5) that was 
present at the start of the recursion. Thus this structure holds at each step. We see further 
that at each step, (Xj^j^ | is obtained from {vj^^^ | by a matrix multiplication involving a 
sum over four spin indices for each choice of and si. At the end, we need a simple matrix 
multiplication to multiply (^sp™)| by \p{{p, f, c}2)spm) ■ From beginning to end, there are 
— 1 steps, where N is the number of final state partons. Thus the amount of computer 
time and storage needed to execute the complete calculation is linear in the number of 
final state partons. One may encounter problems arising from the statistical fiuctuations 
associated with the weight factors, but calculating the weight factors themselves should 
not be a problem. 

4. Properties of the spin decay matrices 

The spin decay matrices ym{s'i, si) are hermitian, 

yg\s[,sir = y^J^{si,4) . (4.1) 
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To see this, we simply note that the starting values of these matrices, <Jsj,<j, , are hermitian 
and that the recursion relation specified by Eqs. (3.11) and (3.11) preserves this property. 
(The coefficients Aik as constructed in Sec. 7 obey Aik{{p}m, s'^, Sm)* = Aik{{p}m, Sm, s'm)-) 
Thus the spin decay matrices have two eigenvectors and have real eigenvalues. The spin 
density matrix for the hard scattering, ({s', s}2|/9({p, /, c}2)spin) in Eq. (2.12) is also her- 
mitian, so the spin weight factor (lspin|/Os^l) = {Y^i^\p{{p, f, c} 2) spin) is real. 

If the recursion relation specified by Eqs. (3.9) and (3.11) were somewhat simpler, 
taking the form 

sJn.Sm s[,Si (4.2) 

si,si)vi 

where C; > 0, then the spin decay matrices ym{s'i, si) would be easily seen to be a positive 
matrix. That is both eigenvalues would be positive and J2 i^'i) ym {s'l, si) u{si) would 
be positive for any vector u. To see this, we simply note that the starting values of these 
matrices, (^sj.s; > are positive. Then we insert eigenvector expansions for the daughter parton 
spin decay matrices, 

y^Hs'i,rsi)=Y,^j,iUii'<)CjArsi) , 

K 

into Eq. (4.2) and form a dot product with u*(s'i) and u{si), where u is an arbitrary vector. 
This gives 

Vm-M' Si) U{si) = ClY^ Xj,l \K,m 
SpSj J,K 

X X] ^(«;)^«({P>/}m,Sm,Sz,Si)0,z(^«)^A',m(^m) 

X X] ^*(«i)^'r(te/}m,S^>'S«,Sz)O,«(s0C/^,m(s',„) • 

(4.4) 



This is a sum of positive eigenvalues Xj^i and XK,m times the square of the absolute value of 
a certain quantity. It is thus positive, showing that y^_i{s'i, si) is a positive matrix. If the 
spin decay matrices are positive, then their product with the spin density matrix for the 
hard scattering, {{s' , s} 2\ p{{p, f,c} 2) spin) in Eq. (2.12), must give a positive spin weight 
factor. 

Since the main term in the recursion relation leads to positive spin decay matrices, one 
may reasonably suspect that the spin weight factor is positive for most events. However, 
since there are additional terms in the recursion relation (3.11), there can be configurations 
in which the spin weight factor has a negative eigenvalue. 
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5. Example of one step in spin evolution 



We can gain some insight into the evolution of the spin decay matrices as we go backward 
from the end of the shower to the beginning. Let us write the spin decay matrix y for the 
mother parton in the form 



(0 



S,,Sl 



(5.1) 



Here is a two component vector normalized to ^ |^;(s)p = 1. It is the eigenvector 
of the matrix y with eigenvalue A-)_ = (1 + P;). The other eigenvector, with eigenvalue 
X- = Yi (1 — Pi), IS orthogonal to ^. We specify which eigenvector is by choosing A+ > A_, 
so that Pi >0 (as long as Y; > 0). It will be helpful to have some names to apply to these 
variables. We can call Yi the spin enhancement factor and Pi the fractional polarization 
of parton I. Then ^i gives the direction of the polarization. This is a modification of the 
conventional language, in which "fractional polarization" and "polarization vector" refer 
to the preparation of a state rather than to its decay pattern, but it seems appropriate 
when we think of the spin structure as evolving backwards, toward the hard interaction. 
For a gluon state, the polarization vector e corresponding to the vector ^ is 



(5.2) 



Here the vectors e^(jo, ,s; Q) arc hclicity eigenvectors, orthogonal to p and to the total final 
state momentum Q. We will be mostly interested in plane polarized states, for which 



C(+l) 



-i<j> 



/V2 



(5.3) 



This makes a real vector, orthogonal to p and Q. It is rotated about p (in the Q = 
frame) through angle (p, starting from a plane that depends on the phase conventions used 
to define e'^(p, s;Q). 

The mother parton spin decay matrix is determined by the spin decay matrices of the 
two daughter partons, which we parameterize as 



y(J^{s'i, si) = Yi \2Pi Ciis'i) ^tisi) + (1 - Pi) Ssi,s, 



ym\^mi ^m) — ^^(Sj„) ^^(s^) + (1 Pm) ^s' ,$„ 



(5.4) 



We consider the example of a final state g ^ g + g splitting. The matrix y\n-i{'^h s'l) 
depends on several parameters. For the sake of illustration, we make arbitrary choices for 
these with the aim of representing a roughly collinear sphtting. With the convention that 
we denote momentum vectors hy p = {p^ , p^ , p'^ , p^) , we take Q = (1000,0,0,0). A gluon 
with momentum along the x axis splits into two gluons with momenta 



= (110.5,110,+10sin(/), -1OCOS0) , 
Pm = {90.Q, 90, -losing, +10 cos 0) . 



(5.5) 
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Here (p gives the azimuthal angle of the sphtting. We will let vary. The splitting includes 
quantum interference in which gluon m can also be emitted by a gluon k with momentum 

Pk = (141.4,140,0,20) . (5.6) 

In this example, the interference contribution is not particularly large because gluon m 
is not soft. We will let the two daughter partons carry 15% linear polarization, in the 
direction = 0.075 vr and = -0.075 vr. Thus Pi = = 0.15. We set Yi = Ym = l. 

We need a definite choice for the function Aik that partitions the interference graphs 
between a term associated with the splitting of parton / and a term associated with the 
splitting of parton k. We adopt the definitions in Eqs. (7.16) and (7.12) from Sec. 7. 

In Fig. 6, we plot the average spin enhancement factor Yi for this splitting as a function 
of the azimuthal angle cj) of the splitting. We see that the enhancement factor varies about 
15% from 1, where 1 corresponds to no enhancement. The mother gluon acquires a net 
polarization in the splitting. The polarization Pi is shown as a function of (p in Fig. 7. The 
mother polarization depends on cj) but its average is about 15%, the input polarization of 
the daughter partons. 

2 . 
1.5 

Yi 1- 

0.5 



— vr TT 

Figure 6: The spin enhancement factor Yi, Eq. (5.1), contributed by a single splitting in which the 
daughter polarizations are 15%. The splitting is roughly coUinear. as described in the text. The 
spin enhancement factor is plotted against the azimuthal angle of the splitting along with a line 
indicating Y; = 1. 

In this example, we used daughter polarizations of 15%. If the daughter polarizations 
are zero, then the enhancement factor, Y;, is 1 for all azimuthal angles (p. How, then, do 
the gluons get polarized? To investigate this, we set the daughter polarizations to 0% and 
plot the resulting polarization Pi of the mother versus the azimuthal angle (p. The result is 
shown in Fig. 8. We see that the mother polarization varies with (p and is generally a little 
larger than 10%. The corresponding polarization vector is approximately in the direction 
of fc_L, the part of pi — pm orthogonal to pi + pm and Q. 

We conclude that at the first step back from the final state, the gluons will become 
polarized. From there on, they will continue to be polarized and their polarizations will 
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0.4 
0.3 

Pi 0-2 
0.1 

— vr 71" 

Figure 7: The mother parton polarization Pi, Eq. (5.1), for tlie spHtting described in the text, 
plotted against the azimuthal angle 4> of the splitting. We also show a line indicating Pi — 0.15, 
the input polarization of the daughters. 

generate spin enhancements in the distributions of the azimuthal angles of the splittings 
beyond what is present in the spin averaged case. Our numerical experiments, like that 
presented here, suggest that the correlations in azimuthal angle between two successive 
splittings are of rather modest size. We expect, however, that correlations among several 
splittings, involving products of the spin factors Yi, can be quite different from 1. 

0.4 
0.3 

Pi 0-2 
0.1 

— TT TT 

Figure 8: The mother parton polarization Pi, Eq. (5.1), for the splitting described in the text but 
with daughter polarizations set to zero. We plot P; against the azimuthal angle (j) of the splitting. 

6. Example of three steps in spin evolution 

We can illustrate the calculation of the spin weight factor for a simple final state shower 
with three splittings. Suppose that shower evolution with spin averaged has produced a 
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shower as illustrated in Fig. 9. The gluon with label has split into gluons with labels 1 
and 2 with the participation of gluon 7 in the soft interference term. Then gluon 1 has 
split into gluons 3 and 4 with the participation of gluon 2 in the interference term. Finally 
gluon 2 has split into gluons 5 and 6 with the participation of gluon 4 in the interference 
term. There are other final state partons that do not participate and are not shown. At 
each splitting, a small amount of momentum is taken from the partons that did not split 
according to the momentum mapping of Ref. [1] . 

Suppose that at the end of the shower the momenta of the final partons, along with 
the momentum of all the final state partons, q, are given in components {E,px,Py,Pz) by 

P3 = (114.564,110,32,1) , 
P4 = (131.852,130,22,-1) , 
P5 = (90.05, 90, 3 cos 0, 3 sin (/>) , 

P6 = (70.0643,70, -3 cos 0,-3 sin 0) , ^^"^^ 
P7 = (131.909,130,-22,-4) , 
q = (1000,0,0,0) . 

Note that the splittings in this example are approximately coUinear and that the azimuthal 
angle of the 5-6 splitting is denoted as (p ^-nd left variable. 

We assume that the spins of the partons at the end of the shower are not measured, so 
that their spin decay matrices are unit matrices. We then compute the spin decay matrices 
of partons 1 and 2. Finally, we use these spin decay matrices to compute^ the spin decay 
matrix y^^^ of parton 0. From this we can use Eq. (5.1) to compute the polarization of 
parton and its spin enhancement factor Yq. If parton came directly from the hard 
interaction, then we would take the trace of y^^'^ with the spin dependent hard matrix 
element to calculate the spin weight factor. Thus the complete spin weight factor contains 
the spin enhancement factor Yq of parton as a factor. 

In Fig. 10, we plot Yq as a function of the azimuthal angle (j) of the 5-6 splitting. We 
see that angles near and vr get higher weights. This illustrates that spin effects act to 
correlate azimuthal angles of splittings even though the spins of the partons at the end of 
the shower are not measured. We also see that the effect with just four active final state 
partons is not large, a little bigger than 10%. We should note, however, that the effects can 
be larger when the depth of the shower is bigger and more final state angles are correlated. 

7. The dipole partitioning function 

Eq. (3.11) contains a function ^;/t({p}m., •^m) that specifies how the two l-k interference 
diagrams are partitioned into separate terms. A fraction Aik is associated with the splitting 
of parton I and comes with the momentum mapping for that splitting, while a fraction A^i 
is associated with the splitting of parton k and comes with the momentum mapping for 
that splitting. We have 

Aik{{p}m,s'^,Sm) + Akl{{p}m,s'„^,Sm) = I ■ (7.1) 
®We again define the function Aik that partitions the interference graphs using Eqs. (7.16) and (7.12). 
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Figure 9: Example of a shower with three sphttings, for which we calculate the spin enhancement 
factor Yq after working back to the first step. 



2 . 



1.5- 




0.5 



— vr 







TV 



Figure 10: The spin enhancement factor Yq for the starting gluon in the mini-shower shown in 
Fig. 9. The spin enhancement factor is plotted against the azimuthal angle (j) of the 5-6 splitting 
along with a line indicating Yq ~ I. 

We can generalize the choice of ^/fc({p}m) given in Ref. [2] in two ways. 
The first choice of Aik is spin independent, 



Aik{{p}m) 



2pi ■ D{p„^,Q) -pk 



Pm -Pl Pm -Pk 



( 



Pl ■ D{pm,Q) -Pl 
iPm-piY 



Ak{{p}m) 



Plk-DiPm,Q)-Plk 
{Pm-Pl Pm-PkY 



) 



(7.2) 



where Q is the total momentum of the final state partons, Pik is the vector 



Pik =Pm-Pi Pk -Pm-Pk Pl 



(7.3) 
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D{pm, QY'^ is the spin sum 



(7.4) 



Pm ■Q (Pm- QY 

and A'li^ is any positive function with 

A'i,{{pU) + A',ii{p}m) = l . (7.5) 

Simple algebra starting from these relations yields the property (7.1) of An^. 

The choice of Aii^. affects the spin-averaged splitting function that appears in Eq. (2.2), 

<^lk{{pJ}m) = CF[Wu-Wik] , (7.6) 

where Wn and Wik are the splitting functions averaged over spin, analyzed in Ref. [2] . As 
in Ref. [2] , it is useful to add and subtract the soft gluon approximation to Wn , 

^eikonal Pl ■ D{pm, Q) ■ Pi 

^" {Pr.-Pl? • ^'-'^ 

The l-k interference function is already constructed using the eikonal approximation that 
gives the soft gluon limit and, with our definition, includes the function Ai^^ 

W,, = A.a.Aj-i;^^h^ . (7.8) 

Pm -Pl Pm- Pk 

Thus we decompose into 

^ik = Cf Wu - WT"-^') + (Wf - Wi,)] . (7.9) 

The term (VI/^'*""^^ — Wik) includes the soft singularity and the softxcollinear singularity, 
while (Wu — Wu ) has only a coUinear singularity. After a little bit of algebra, one 
obtains 

Y^eikonal Yf7 a A' ^Ik ' D (pm , Q) ■ Plk . ./ ~^lk ,"7 i n\ 

[Pm-piPm-Pk) [Pm-piPm-Pk) 

The last equality here follows from the fact that pm ■ Pik = 0. This has the feature that 
~ Wik) is positive as long as A[f^ is positive. This is important for constructing 
the spin- averaged shower as a Markov process. 
The choice of in Ref. [2] was 

./ /jn ^ [Pm-Pk? Pl ■ D{Pm,Q) - Pl 

AkiiPim) = — — — — — — — 7—- . (7.11) 

[Pm ■ PkY Pl ■ D{p.m, Q) ■Pl + [Pm ■ PlY Pk ■ D[p„i, Q) ■ Pk 

With this choice. An. = A'j^^. This choice has three good features. First, it vanishes 
Pm -Pfc — > in such a way that there is no singularity in {wf^^"^^"^ — Wik) when pm becomes 



- 21 - 



coUinear to pk- Second, it is invariant under rescaling of pm- Third, it is also invariant 
under rescaling of pi and of pk so that, for the case of massless partons I and fc, it is a 
function only of the angles of the partons. 
Another possible choice is 



This choice has the same three good properties, but it is simpler. A third possible choice 
is 



This has two of the good properties but is not invariant under rescaling of pi and pk, so 
that the partitioning depends on the energies as well as the angles of partons / and k. This 
choice does have the advantage of being the simplest. It is analogous to the partitioning 
factor in the Catani-Seymour scheme for dipole subtractions [9]. 

One can also choose a spin dependent partitioning factor. This is particularly useful 
in the case, investigated in this paper, that we take the leading color approximation. We 
rewrite Eq. (3.11) by adding and subtracting the eikonal approximation for the direct l-l 
splitting graph and by using a spin dependent Aik as 






(7.13) 





(7.14) 



X V2,l{{P^ hm, Sm, Si, Si) V^idp, f]m, s'^, s'l, s'{) 
- Vz,l{{P, f}m, Sm, Si, Si) V^i{{p, f}m, s'm, s'l, S;) 




Aik{{p}m,Sm,S, 



e{Sm)*-pl e{Kn)-Pk + £{Sm)*-pk £{s'm)-Pl 



'm 



Pm.-Pl Pm-Pk 



Here we have adopted the shorthand notation 





(7.15) 



^The singularity wlien pm becomes collinear to pi is assigned to {W] 



rcikonal 
kk 



Wkl). 
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One can set 

Pm-Pl Pm-Pk 



^lk{{p}m, Sr^ 



e{sm)*-pi e{s'^)-pk + e{sm)*-Pk e{s'^)-pi 



{Pm-piY " {Pm-PlPm-PkY 



(7.16) 



The matrix A'^^, should be positive and obey A'^^ + yl^^ = 1 . Thus one can take one of the 
choices given above for it. Then 

Aik{{p}m,s'^,Sm) + Akl{{p}m,S^,Sm) = 1 ■ (7.17) 

With this form for A/^, the last term in Eq. (7.14) simplifies so that we obtain 
v'^^ is', sA - '^F'^'^i-^^"') V 'u('")fs' s ^ Vu«fs', sA 

ym-l\^h^l)— , s tlm V*m'*mJ / , Um V^h 

'^Ikiipjjm) J, s 5' 5, 

X jv/ ■5/,S«)^r({P, /}m,Sm>S/)'Si) 

. r e £(^m)*-Pf £(g'm,)-P^ 

+ e{le{l,...,m-l}Ji = fm = g) 

V2,l{{P^ f}m, Sm, k, Si) vl^i{{p, f}^, s'm, s'l, s'l) 
- ^3,/({P, f}m, Sm, h, Si) /({p, f}m, S;, s'l) 

The last term in Eq. (7.18), constructed from the eikonal approximations for the direct and 
interference graphs, contains the soft and the soft x coUinear singularities. The rest of the 
contributions together represents the pure coUinear singularities. In Eq. (7.18), the physics 
of the soft and soft x coUinear term is clear. The polarization of parton I is unchanged and 
the plane of the decay (as specified by Pik) is preferentially aligned with the polarization 
of the emitted soft gluon. 

We also need the spin averaged splitting function ^ik with this choice of Aik- Rather 
than start with its definition, we can use Eq. (7.14). According to Eq. (2.5), if we insert 
unit matrices for y^\s'^, Sm) and ym{s\, si), the output y^^}_i{s[, si) must have trace equal 
to two. This fixes the factor in Eq. (7.14). This calculation gives Eqs. (7.9) and (7.10) 
for with the factor that was used for the spin dependent Aik-, Eq. (7.16). 

It is of interest to examine the soft gluon spin factor from Eq. (7.18), 



^Ik 

This is normalized to ^sF{s, s) = 1, so that it represents the contribution to the spin 
weight function from a single splitting in which the new parton m is a very soft gluon. For 
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this purpose, we assume that partons I and k are parts of a narrow jet with momentum 
in approximately the z-direction in the rest frame of Q, the total momentum of all the 
final state partons. Then we can parameterize pi, pk and Pm by using two dimensional 
transverse vectors 0;, 0fc and 6m as, with the notation p = (j)^,p,p^), 



Pi = Ei{i,ei,^i-ef) , 

Pk = Ekil,ek,^Jl-el) , (7-20) 

Pm = Em (1, Om, \/ 1 — ^m) ■ 

We will assume that \6i\ <C 1 and \6k\ ^ 1. Then, because of destructive interference, 
the important integration region for 9m is \0m\ ^ 1. We can then make small angle 
approximations everywhere. For the polarization vectors we can approximate 

e^{0,s,-em-e) (7.21) 

with |e| ~ 1. Then pm • e ~ 0. With these approximations 

[jOl - Om? {Ok - Om) ■ ejpm, Sm] Q)* - {Ok - Ofnf {Ol - dm) ' sjpm, s'm', QT (^"^^^ 
{Oi - OmYiOk - - 0kY 

We can combine F{s\ s) with the spin decay matrix for parton m in the form of 
Eq. (5.1) 

Vm^ i.^m^^m) = Ym\^Pm^m{Sm) ^m{^m) ~\~ ~ Pm) Ss'^^Sm ■ (7.23) 

In Fig. 11, we display 

{F)= E y'^:'^{SmrSm)F{s'mrSm) (7.24) 

as a function of the angles {9m,x,Gm,y) of the soft gluon m. We have taken Oi = (0.1,0) 
and Ok = (—0.1, 0). We have set Ym = 1 with polarization Pm = 0.15 in the direction 

^L{s'm)<Pm.,s'm;Q) = (1,0) . (7.25) 

s 

We see that for this polarization choice, the spin factor gives positive interference for Om. 
in the region between Oi and 0^. This is in addition to the positive interference in the 
spin-averaged cross section in this region that we took note of in Ref. [2] . For polarization 
in the (0, 1) direction, one has the opposite effect, (F)(o,i) = 1 — (-F)(i q)- 



8. Conclusions 

In Ref. [2], we have seen how to formulate a leading color, spin averaged shower version of 
the general parton shower formalism of Ref. [1] in a fashion that could be implemented in 
a standard style of calculation as a Markov process. However, this approximation leaves 
out potentially important effects related to spin and color. In this paper, we have seen how 
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Figure 11: The soft gluon spin weight function (F) from Eq. (7.24), plotted against the emission 
angle {Om.x,0m.,y) of the soft gluon. Parton I is at 6i = (0.1,0) and parton k is at 0^ = (—0.1,0). 
The soft gluon is assumed to be 15% polarized in the (1,0) direction. 

one might add back the effect on the parton momentum distributions of the parton spin 
correlations, while still working in the leading color approximation. 

We assume that the spins of the partons that appear at the end of the shower are not 
measured by the measurement function applied to the state produced after hadronization 
of these partons. Then the probability to get a given shower history and set of parton 
momenta at the end of the shower is the same probability as for a spin averaged shower 
but times a certain spin weight function. We have seen that this spin weight function can 
be obtained from a straightforward computation that takes an amount of computer time 
that is linear in the number of partons. 

If a parton shower event generator produces events with weights Wi instead of events 
with weight 1, there is a potential to degrade the numerical convergence rate of the cal- 
culation of a desired cross section: it may take more Monte Carlo events to produce a 
result with the same statistical accuracy. This can happen if the weights have a large vari- 
ance. For instance, suppose that for a certain observable of interest the weights average 
to W = 0.6, indicating that the expectation value of this observable is only 60% as large 
as it would be for a spin- averaged shower. If almost all of the weights lie between —Ixw 
and 3 X w, there is no real loss of accuracy, but if the weights are typically spread over the 
range —100 x w < Wi < +100 x w, then many more Monte Carlo events will be needed to 
produce the same accuracy as with all equal weights. In this event, the simplest approach 
would be to calculate the spin weight factor corresponding to just the first N splittings 
after the hard interaction. The true answer is then found by taking the limit — > oo. 
The maximum value of N that could be used would be limited by the computer power 
available. It can well be that, for a final state observable of interest, iV = 10 is enough. 

On the other hand, if some final state observables are found to be very sensitive to spin, 
then, for those observables, larger values of N than are practical might needed. In that 
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case, one could modify the algorithm for generating events by choosing azimuthal angles 
with a probability p' = F x pQ, where po is the probability for generating events in the 
spin averaged case and -F is a new factor. For example, one could follow Pythia and make 
use of the correlations given in Ref. [6] between the azimuthal angle of a parton splitting 
and the azimuthal angle of the splitting that produced the mother parton. Whatever 
method is used, one wants F > 1 for the kinds of events for which the spin weight factor 

= (ispin |pg in) is larger than 1 and F < 1 for the kinds of events for which u;o < 1. 
This would give a new spin weight factor w' = wq/F, so that the new weight factor more 
nearly approximates 1. It is not necessary to get F to exactly match wq. One simply has 
to arrange that wq/F is never very large compared to its average value. 

We have given formulas for the calculation of the spin weight function using the split- 
ting functions that appear in the general formulation of Ref. [1] . However, the basic method 
is of quite general applicability. To include spin in a parton shower Monte Carlo based on 
averaging over spins, one needs spin dependent splitting functions with two properties: 1) 
their soft and collinear limits match QCD matrix elements and 2) their spin averages are 
the splitting functions used in the existing parton shower Monte Carlo. For instance, one 
could define spin dependent dipole splitting functions that generalize those of Catani and 
Seymour [9] for use in a dipole shower [10]. In addition, it should be straightforward to 
define spin dependent antenna splitting functions [11] for use in an antenna shower [12]. 
With spin dependent splitting functions, one can generate the spin weight functions by 
working backwards from the final state, as we have seen in this paper. 
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